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RESEARCH  IN  NONLINEAR  PARTIAL  DIFFERENTIAL  EQUATIONS 

AND  BIFUCATION  THEORY 
Grant  #  AFOSR-86-0218 
FINAL  REPORT  -  December  22,  1987 

Dr.  David  H.  Wagner 
University  of  Houston 
Department  of  Mathematics 
4800  Calhoun 
Houston,  TX  77004 

Significant  progress  has  been  made  under  the  sponsorship  of  AFOSR  Grant  number 
AFOSR-86-0218,  on  the  mathematical  theory  of  combustion  waves  at  the  modelling  level  of  the 
Navier  Stokes  equations  for  a  reacting  gas.  In  particular,  simple,  "clean,"  and  relatively  sharp 
criteria  have  been  found  which  ensure  the  existence  or  non-existence  of  strong  plane  detonation 
wave  solutions  to  the  reacting  compressible  Navier  Stokes  equations,  under  the  assumption  of  a 
simple  one  step  reaction  with  generalized  Arrhenius  type  kinetics  and  arbitrary  Lewis  numbers. 

These  results  give  a  mathematical  proof  of  the  structural  stability  of  ZND  detonation  waves 
under  the  influence  of  heat  conductivity,  viscosity,  and  chemical  diffusion.  It  has  also  been 
possible  to  identify  parameter  regimes  in  which  variables  such  as  the  temperature  of  the  gas  are  or 
are  not  strictly  increasing  throughout  the  wave.  Rigorous  verification  has  been  given  of  a 
phenomenon  discovered  in  machine  calculations  by  Majda,  Colella,  and  Roytburd,  namely  that 
when  the  product,  of  the  heat  conductivity  and  the  strong  detonation  burned  state  value  of  the 
Arrhenius  function,  is  large,  then  the  profile  of  the  strong  detonation  wave  breaks  up  into  two 
separate  waves.  These  waves  look  like  a  weak  detonation  wave  followed  by  a  non-reacting  shock 
wave.  Although  in  most  applications  these  waves  are  considered  to  be  non-physical,  they  present 
a  real  problem  for  machine  calculations.  The  reason  for  this  is  that  a  small  amount  of  artificial 
diffusion,  in  combination  with  a  physically  very  large  burned-state  value  for  the  Arrhenius  function 
can  easily  trigger  these  misleading  waves. 

The  details  of  these  results  arc  described  in  the  attached  preprint: 

"The  existence  and  behavior  of  viscous  structure  for  plane  detonation  waves,”  David  H. 

Wagner,  submitted  to  SI.VM  J  Math.  Anal. 

In  addition,  the  paper  i.iied  below  has  been  published,  which  shows  the  right  way  to 
rigorously  make  the  transfomiaiion  fn»m  Fulerian  coordinates  to  Lagrangian  coordinates  for  weak, 
shock-wave  solutions  to  the  equaiions  of  mviscid  gas  dynamics.  This  result  will  certainly  become 
a  part  of  many  graduate  students  training  It  also  gives  some  valuable  insights  into  the  correct 
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Tcprescntation  of  &  vacuum  in  Lagrangian  coordinates.  Reprints  have  been  submitted  in  an  earlier 
report. 


"Equivalence  of  the  Euler  and  Lagrangian  equations  of  gas  dynamics  for  weak  solutions, " 
David  H.  Wagner,  Jour.  Diff.  Eq.,  1987,  118-136. 


THE  EXISTENCE  AND  BEHAVIOR  OF  VISCOUS  STRUCTURE  FOR 

PLANE  DETONATION  WAVES 


David  H.  Wagner* 

Department  of  Mathemadcs 
University  of  Houston 
Houston,  TX  77004 

ABSTRACT 

.  We  prove  a  necessary  condition  and  a  sufficient  condition  for  the  existence 
of  steady  plane  wave  solutions  to  the  Navier  Stokes  equations  for  a  reacting  gas. 
These  solutions  represent  plane  detonation  waves,  and  converge  to  ZND  detonation 
waves  as  the  viscosity,  heat  conductivity,  and  species  diffusion  rates  tend  to  zero. 
We  assume  that  the  Ptandti  number  is  3/4 ,  but  we  permit  arbitrary  Lewis  numbers. 
We  make  no  assumption  concerning  the  activation  energy. 

We  show  that  the  stagnation  enthalpy  and  the  entropy  flux  are  always 
monotone  for  such  solutions,  and  that  the  mass  density  and  pressure  are  nearly 
always  not  monotone,  as  predicted  by  the  ZND  theory. 

In  certain  parameter  ranges,  typically  that  of  large  diffusion,  many  of  these 
waves  have  the  appearance  of  a  weak  detonation  followed  by  an  inert  shock  wave. 
This  confirms  a  phenomenon  observed  in  numerical  calculations  and  in  a  model 
system  by  Colella,  Majda,  and  Roytburd. 

1.  INTRODUCTION 


Detonation  waves  are  compressive,  exothermically  reacting  shock  waves.  One  of  the 
curiosities  of  combustion  theory  is  that  there  also  exist  expansive  "shock  waves"  known  as 
deflagration  waves,  which  will  not  be  discussed  in  this  paper.  We  will  give  a  mathematically 
rigorous,  but  simple,  discussion  of  the  viscous  structure  of  plane  detonation  waves. 

We  begin  with  a  brief  discussion  of  the  inviscid  theory,  known  as  the  Chapman  Jouguet 
(CJ)  theory.  If  we  assume  that  the  thickness  of  the  reaction  zone  is  zero,  if  we  neglect  all  diffusion 
effects  such  as  viscosity,  heat  conduction,  and  diffusion  of  species,  and  any  external  forces  such 
as  gravity,  and  if  we  look  for  steady  plane  waves,  then  we  obtain  the  following  system  of 
differential  equations: 

(1.1)  (a)  (pu);^  =  0, 

(b)  [pu2  +  p(p,T)],  =  0, 

(c)  [[p(u2/2  +  e(p,  T,  Y))  +  p(p,  T)]u];^  =  0, 

(d)  (puY)x=  -  puY- 5(x-xo). 


Here  x  is  a  space  coordinate  in  the  direction  normal  to  the  wave,  xq  is  the  location  of  the  wave, 
and  p,  T,  u,  p,  e,  and  Y  are  the  mass  density,  temperature,  x-component  of  velocity,  pressure, 
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specific  internal  energy,  and  mass  fraction  of  the  reactant,  respectively.  As  is  standard  practice,  we 
have  represented  the  extremely  complicated  chemical  reaction  by  a  simplified,  one-step  chemistry: 
reactant  — >  product.  From  (1.1a)  we  see  that  the  mass  flux,  pu,  has  a  constant  value;  we  denote 
this  value  by  m.  The  fluxes  of  momentum,  (1.1b),  and  energy,  (1.1c)  are  also  constant;  from  this 
fact  we  obtain  the  Rankine-Hugoniot  conditions  for  a  shock  wave,  which  in  the  inviscid  theory  is 
represented  by  a  jump  discontinuity  in  the  unknowns.  The  difference  between  inert  gas  dynamics, 
and  ihe  exothermic  reactive  theory  discussed  here,  lies  in  the  fact  that  Y  varies  from  a  positive 
value  on  the  unbumed  side  of  the  wave,  which  we  take  to  lie  on  the  right  side,  to  a  zero  value  on 
the  burned,  or  left  side.  Because  the  internal  energy  e  depends  on  Y,  the  change  in  Y  causes 
the  classical  Hugoniot  curve  (the  solution  locus  of  (1.1c))  of  gas  dynamics  to  move.  As  a 
consequence,  we  find  that,  for  a  given  value  of  m,  a  given  shock  state  on  the  left  may  now  be 
coimected  by  a  shock  wave  to  two  possible  states  on  the  right,  except  for  certain  critical  values  of 
m  for  which  there  is  a  uruque  burned  state  —  the  Chapman  Jouguet  point  —  see  Bg.  1.  In 
P 


addition,  the  curve  of  possible  burned  states,  parameterized  by  m,  has  two  components.  One 
component,  corresponding  to  compressive  waves,  is  called  the  detonation  branch,  and  the  other 
component,  corresponding  to  expansive  waves,  is  called  the  deflagration  branch.  By  way  of 
contrast,  in  an  inert  gas,  for  a  given  value  of  m,  a  state  is  usually  connected  to  only  one  state  on  the 
right,  and  the  curve  of  possible  terminal  shock  states  is  usually  connected. 
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The  combustive  shock  waves  of  the  CJ  theory  are  classified  as  follows.  A  wave 
connecting  the  onbumed  state  to  the  closer  detonation  point  is  called  a  weak  detonation  wave,  and 
a  connection  to  the  farther  detonation  point  is  called  a  strong  detonation.  A  detonation  wave 
terminating  at  the  Chapman  Jouguet  point  is  called  a  Chapman  Jouguet  detonation.  Deflagration 
waves  are  similarly  classified.  Weak  detonations  are  generally  not  observed  for  the  exothermic, 
irreversible  reactions  considered  here,  but  may  be  observed  in  other  contexts  [FD].  Strong 
deflagrations  violate  the  second  law  of  thermodynamics  and  are  unphysical. 

The  CJ  theory  for  detonation  waves  is  useful  for  deriving  the  Rankine  Hugoniot 
conditions,  and  for  classifying  the  types  of  wave.  However,  this  theory  is  physically  flawed, 
because  in  reality  the  reaction  zone  is  much  thicker  than  the  shock  layer.  This  is  due  to  the  fact  that 
the  chemical  reaction  depends  on  molecular  collisions,  and  requires  a  distance  much  longer  than  the 
mean  free  path  to  achieve  significant  completion.  The  shock  layer,  however,  has  been 
experimentally  observed  to  be  several  mean  free  paths  thick.  Consequently  the  appropriate  inviscid 
model  is  the  one  developed  independently  by  Zel’dovich,  von  Neumann,  and  D6ring[Z,  Nl,  N2, 
D],  and  which  is  known  as  the  ZND  model.  In  this  model  equation  (l.ld)  is  replaced  by  a  similar 
equation,  but  with  a  finite  reaction  rate: 

(l.ld*)  (puY\  =  -  r(p,  Y,  T) 

For  our  purposes  it  is  reasonable  to  assume  that  the  reaction  rate  function  r  is  continuous, 
non-negative,  and  monotone  in  each  variable.  Our  mathematical  treatment  will  require  that  we 
assume  that  r  vanishes  whenever  the  temperature  T  is  less  than  a  given  ignition  temperature  Tj. 
For  a  known  reaction  rate  r  one  can  solve  (l.la,  b,  c,  d')  explicitly;  the  only  detonation  wave 
solutions  are  strong  or  CJ  detonations.  These  waves,  which  are  known  as  ZND  waves,  begin 
with  a  jump  discontinuity  which  is  an  inert  shock  wave.  This  shock  wave  heats  the  gas  above  the 
ignition  temperature;  the  reaction  proceeds,  with  the  velocity  and  temperature  following  a  curve  of 
equilibrium  states  for  (1.1b,  c),  parameterized  by  Y.  One  of  the  interesting  features  of  these 
waves  is  the  peak  in  the  pressure  and  density  which  is  known  as  the  von  Neumann  spike.  By  way 
of  contrast,  in  inert  shock  waves  the  pressure  and  density  usually  monotone  [Gi]. 

The  equations  of  inert,  inviscid.  non-heat-conducting  gas  dynamics  are  an  example  of  a 
nonlinear  hyperbolic  system  of  conservation  laws.  In  the  theory  for  such  systems  it  is  standard 
practice  to  set  admissibility  critena  to  distinguish  physical  from  unphysical  shock  waves.  One  of 
the  criteria  in  which  much  faith  is  put  is  to  accept  a  shock  wave  as  physical  if  it  is  structurally 
stable.  A  shock  wave  is  structurally  stable  if  it  is  the  limit  of  solutions  to  models  which  include 
more  physical  effects,  such  as  viscosity  and  heat  conduction,  as  these  model j  tend  to  the  original 
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inviscid  model  in  which  these  effects  are  neglected.  For  steady  plane  detonation  waves  one  may 
consider  the  effects  of  viscosity,  heat  conduction,  and  species  diffusion,  to  obtain  the  (steady) 
reacting  compressible  Navier  Stokes  equations: 


(1.2)  (a)  (pu)^  =  0. 

(b)  [pu2  +  p(p,  Dljt  =  (pux)x, 

(c)  [[p(u2/2  +  e(p.  T,  Y))  +  p(p,  T)]u], 

(d)  (puY);,=  (pDYA-r(p.T,Y). 


=  (XT,),  +  (puu,),  +  (qpDY,), , 


Here  p,  is  the  coefficient  of  bulk  viscosity,  X  is  the  heat  conductivity,  D  is  the  diffusion  rate  for 
the  reactant,  and  q  is  the  difference  in  the  heats  of  formation  of  the  reactant  and  the  product  [Wi]. 

In  this  paper  we  prove  a  necessary  condition  and  a  sufficient  condition  for  the  existence  of 
heteroclinic  solutions  of  (1.2)  which  extend  from  an  unbumed  state  at  x  =  -<»  to  the  strong 
detonation  point  at  x  =  <».  These  conditions  also  apply  to  the  Chapman  Jouguet  detonation.  See 
(4.8)  and  (5.3).  For  simplicity  we  restrict  our  attention  to  the  case  where  the  Prandtl  number  is 
3/4  (p,  =  X/Cp).  In  the  limit  as  X,  p,  and  D  tend  to  zero  (with  other  parameters  fixed)  these 
solutions  tend  to  the  ZND  wave.  Thus  the  ZND  wave  is  structurally  stable  to  this  particular 
perturbation  of  the  model. 

For  all  of  these  solutions  the  stagnation  enthalpy  H  =  CpT  +  u^/2  is  monotone,  as  is  the 
entropy  flux:  mS  -  XT,/T.  For  most  of  the  strong  detonation  waves  the  density  and  the  pressure 
attain  their  maxima  in  the  interior  of  the  wave;  this  corresponds  to  the  von  Neumann  spike  which 
occurs  in  the  ZND  wave.  However  for  a  cenain  parameter  range,  namely  whenever 

J  y 

where  L  is  the  Lewis  number,  M*  and  T*  are  the  Mach  number  and  tempjerature  at  the  strong 
detonation  point,  and  k(()(T*)  is  the  reaction  rate  as  defined  in  (2.9),  there  exists  a  continuum  of 
solutions  which  look  like  a  weak  detonation  followed  by  a  gas  dynamic  shock  wave.  For  these 
waves  the  pressure  and  temperature  are  monotone.  This  pathological  behavior  has  been  noted 
before  in  [Wo2,  FD,  LL,  HoSt|,  and  m  numerical  computations  of  solutions  of  the  time  dependent 
Navier  Stokes  equations  for  a  reacting  gas  (CMR).  In  these  numerical  computations  it  was 
observed  that  the  weak  detonation  shix  k  wave  solutions  are  dynamically  stable  as  solutions  of  the 
time  dependent  equations  in  one  ^p.lwc  dimension. 

Since  Zel'dovich,  von  Neumann,  and  Doring  [Z,  Nl,  N2,  D]  described  the  typical  plane 
detonation  heuristically  as  an  men  Nhixk  wave  followed  by  a  deflagration,  there  have  been  a 


number  of  papers  on  the  structure  problem  for  plane  detonations.  A  common  assumption  has  been 
that  the  Prandtl- number  is  3/4  and  that  the  Lewis  number  is  1.  Under  these  assumptions 
Hirschfelder  and  Curtiss  gave  a  good  analysis  of  the  behavior  of  structure  profiles  [HC],  and 
Wood  proved  the  existence  of  structure  for  "small"  reaction  rates  [Wol].  The  approach  taken  here 
has  much  in  common  with  Wood’s,  except  that  we  have  been  more  precise  in  our  analysis,  our 
results  are  stronger  and  more  general,  and  we  give  explicit  and  fairly  sharp  statements  of  just  how 
small  the  rate  parameter  must  be. 

A  typical  expression  for  the  reaction  rate  is  given  by  the  Arrhenius  law: 


r(p,T,Y)  =  kpYe* 


where  0  is  the  activation  energy  E/R.  In  mathematical  combustion  theory  it  is  standard  practice 
to  simplify  models  such  as  (2)  by  taking  an  appropriate  distinguished  limit  as  k  and  0  tend  to 
infinity.  This  has  the  effect  of  reducing  the  thickness  of  the  reaction  zone  to  zero.  Bush  and 
Fendell  [BF]  gave  a  description  of  CJ  detonations  using  asymptotic  expansions  in  the  limit  of 
infinite  activation  energy.  Stewart  and  Holmes  proved  the  existence  of  viscous  structure  for  (2), 
assuming  a  large,  finite  activation  energy  [HoSt].  Lu  and  Ludford  gave  a  simplified  analysis  of 
weak,  strong,  and  CJ  detonations  in  the  infinite  activation  energy  limit  [LL]. 

It  is  desirable  to  understand  the  existence  of  structure  for  plane  detonation  waves  with  no 
assumptions  concerning  activation  energy,  independent  of  any  desire  for  marhematical  generality. 
Because  (under  the  ZND  hypothesis)  the  reaction  zone  is  much  thicker  than  the  shock  layer,  one 
should  study  this  problem  for  activation  energies  which  are  small  compared  to  the  reciprocals  of 
the  heat  conductivity,  viscosity,  and  species  diffusion  rate. 

Gardner  proved  the  existence  of  travelling  plane  detonation  wave  solutions  to  the 
Lagrangian  reacting  compressible  Navier  Stokes  equations  [Ga].  He  made  no  assumption  on  the 
Lewis  and  Prandtl  numbers.  However,  he  omitted  the  species  diffusion  term  from  the  energy 
balance  equation,  and  this  term  is  not  usually  neglected.  It  may  be  that  the  effect  of  this  term  on  the 
solution  is  very  small.  However,  we  will  demonstrate  in  this  paper  that  inclusion  of  this  term 
permits  a  much  more  natural  treatment  of  the  problem,  and  better  bount  s  on  the  solution. 

The  remainder  of  the  paper  is  organized  as  follows.  In  §  2  we  explain  our  assumptions  in 
more  detail,  and  we  show  that  under  these  assumptions,  and  the  assumption  that  the  Lewis  number 
is  one,  that  (1.2)  is  reduced  to  a  system  of  three  differential  equations.  This  material  has  been 
extracted  and  specialized  from  [Wi]  It  is  included  for  clarity  and  completeness  of  exposition.  In 
§  3  we  present  the  simple  topological  argument  that  proves  the  existence  of  viscous  structure  for 
steady  plane  detonation  waves.  In  §  4  we  prove  the  estimate  that  makes  the  topological  argument 
work,  and  we  conclude  this  section  with  a  sufficient  condition  for  the  existence  of  structure.  In  §  5 
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wc  prove,  using  simple  energy  estimates  on  the  stagnation  enthalpy,  a  necessary  condition  for  the 
existence  of  structure;  if  this  condition  is  not  satisfied  then  viscous  structure  does  not  exist.  In  §  6 
we  give  the  generalization  to  arbitrary  Lewis  numbers.  In  §  7  wc  discuss  the  behavior  of  the 
solutions  for  various  parameter  values.  In  §  8  we  present  a  rigorous  discussion  of  the  ZND  limit. 
We  conclude,  in  §  9,  with  a  proof  that  the  entropy  flux  is  monotone. 

2.  Reduction  to  Three  Equations 

We  make  the  basic  thermodynamic  assumption  that  both  the  reactant  and  product  satisfy  the 
same  ideal  gas  law,  and  differ  only  in  their  heats  of  formation.  Thus  the  pressure  is  independent  of 
Y: 

(2.1)  p  =  RpT. 

We  further  assume  that  the  internal  energy  depends  linearly  on  Y,  so  that  we  have: 

(2.2)  e  =  CvT+qY, 

where  Cy  is  the  specific  heat  at  constant  volume.  These  assumptions  are  probably  not  essential  to 
the  results  that  follow,  and  one  could  probably  replace  them  with  more  qualitative  conditions 
similar  to  Weyl's  conditions  for  the  equation  of  state  for  an  inert  gas.  However,  these  assumptions 
are  essential  to  the  simplifications  that  follow. 

Observe  that  (1.2)  can  be  integrated  once  to  yield  pu  =  m  =  constant.  If  for  any  unknown 
U  we  let  U+  be  the  limit  of  U  as  x  tends  to  ±  oo,  we  have: 

(2.3)  a.  puj^  =  m(u-u^  +  p-p^ 

b.  XT  +  puu^  +  qpDY^  =  m[c^(T  -  T^  )  +  q(Y  -  Y^  )  +  ^u^  -  u^^)  +  R(T  -  T^)] 

c.  (pDY^)^  =  mY^  +  r(p,  Y.  T) 

Let  H  =  CpT  +  u^/2  be  the  stagnation  enthalpy.  Here  Cp  =  Cy  +  R  and  is  the  specific  heat 
at  constant  pressure.  Then  (2.3b)  may  be  rewritten  as: 

(2.4)  (Vcp)Hx  +  (|i  -  Xycp)(u-/2)x  =  m[H  -  H+ +  q(E  -  e±)]  , 

where  e  =  Y  -  pDY^/m  is  a  reaction  progress  variable.  When  the  Prandtl  number  is  3/4,  then 
p.  =  X/Cp  and  wc  have: 


rt* 

*•( 

m 


K: 


(2.5)  =  H-H^  +  q(e-e,). 
p  - 

It  is  convenient  to  let  y  satisfy 
me 

(2.6)  ^  =  — 

dx  X 

so  that  the  left  hand  side  of  (2.5)  becomes  H^.  The  system  (2.3)  reduces  to  a  system  of  four  first 
order  equations: 


(2.7)  a.  mu  =m(u-u^  +  mR[—  — 


■u  u_. 


b.  H  =H-H^  +  q(e-e^) 


c.  e,  =  --^r(p,Y,T) 
m  c 

p 


d.  =  _^(Y-e) 

^  pDCp 

Note  that 

pDc 

(2.8)  (H  -  H^)y  + - £  qYy  =  H  -  +  qY. 

Thus  if  the  Lewis  number,  L  =  XJpDCp  ,  is  one,  we  see  that  the  quantity  H  -  H.  +  q(Y  -  Y.) 
satisfies  the  differential  equation  f '  =  f  and  can  be  bounded  only  if  it  is  identically  zero.  As  we  arc 
interested  only  in  bounded  solutions,  we  may,  in  this  case,  restrict  our  attention  to  the  plane 
H  -  H  +  q(Y  -  Y_)  =  0,  and  (2.7)  reduces  to  a  system  of  three  equations:  (2.7a,  b,  c)  with  Y 
replaced  by  (H+  -  H)/q.  We  will  not  actually  need  to  assume  that  the  Lewis  number  is  1, 
however,  this  case  is  easier  to  understand. 

Some  of  the  results  concerning  this  system  arc  easier  to  interpret  if  we  replace  u  by  the 
specific  volume  V  =  1/p  =  u/m.  In  this  case  (2.7a)  becomes: 


(2.7a') 


V.  -  V- V++- 


R  (T  T, 


nv  T7 


With  the  exception  of  our  discussion  of  entropy,  we  will  restrict  our  attention  to  the  system 
(2.7a’,  b,  c)  =  (2.7)  for  the  remainder  of  this  paper. 
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Following  standard  practice,  we  will  resolve  the  cold  boundary  difficulty  by  means  of  an 
ignition  temperature  assumption.  The  cold  boundary  difficulty  consists  of  the  fact  that  the 
Arrhenius  reaction  rate  (1.3)  does  not  vanish,  but  is  merely  very  small,  at  the  unbumed  state. 
Thus  no  solution  of  (2.8)  can  tend  to  the  unbumed  state  as  x  — >  ».  Clearly  this  problem  stems 
more  from  oiu’  unphysical,  infinite  domain  than  from  any  flaw  in  the  reaction  rate.  However,  as  is 
customary,  we  will  resolve  this  problem  by  modifying  r  so  that  it  is  zero  for  T  <  T, ,  where  Tj 
is  the  ignition  temperature,  and  is  chosen  to  be  greater  than  the  unbumed  temperature  T. . 

Accordingly  we  will  consider  reaction  rates  of  the  form: 

(2.9)  r(p,  Y,  T)  =  kpY  (j)(T) , 

where  (j>  is  a  Lipschitz  monotone  function  of  T,  which  vanishes  for  T  <  Tj . 

In  the  next  three  sections  we  will  usually  assume  that  the  Lewis  number  is  one,  for  ease  of 
understanding.  In  §  6  we  will  explain  how  to  generalize  to  arbitrary  Lewis  numbers. 

3.  A  picture  is  worth... 

We  now  describe  a  region  in  (V,  H,  e)  space.  The  topological  properties  of  the  flow  for 
(2.7)  in  this  region  will  imply  the  existence  of  the  desired  structure  profiles,  and  will  yield  other 
properties  as  well. 

For  any  bounded  solution  of  (2.7),  on  which  some  quantity  is  strictly  monotone,  (V,  H,  e) 
must  tend  to  a  rest  state  of  (2.7)  as  y  tends  to  ±  «x>.  From  (2.7c)  we  see  that  r(p,  Y,  T)  must 
vanish  at  any  rest  state.  For  our  modified  kinetics  this  can  only  happen  if  Y  =  0  or  T  <  Tj.  Rest 
states  satisfying  Y  =  0  are  possible  burned  states;  those  satisfying  T  <  Tj  are  possible  unburned, 
or  fresh  gas  states.  We  choose  to  place  the  burned  state  at  y  =  +<»  and  the  unbumed  state  at 
y  =  -oo.  Thus  Y+  =  0  and  e+  =  (Y  -  pDY;^/m)+  =  0.  Note  that  with  this  orientation  it  is  natural 
to  expect  that  m  =  pu  is  positive. 

At  the  unbumed  state,  Y  has  a  given  value,  Y. ,  which,  since  Y  is  a  mass  fraction,  is 
naturally  bounded  by  one.  Since  E.  =  Y.,  we  have,  by  (2.7b),  R,.  -  H_  =  qY. ,  which  is  the  total 
heat  per  unit  mass  released  by  the  reaction.  Then  (2.7a’)  yields 


(3.1)  = 


R  T^  T 


=  (V^-V)  1  - 


VJ 


Rl  R  (H  +  qY.  H. 


Let  y  be  the  ratio  of  specific  heats.  Cy/Cp  =  (Cp  -  R)/Cp .  Then  (3.1)  becomes: 


I 


(3.2) 


V2- 


V  , 

m^(Y+-l)V. 


lv++  (h.  -h  qY.)  =0. 


If  we  solve  this  for  V^.  we  obtain 


(3.3) 


2(Y-l)qY. 

(y<-l)u? 


where  =  u^/c^  =  u2/(7p/p)  is  the  square  of  the  Mach  number,  and  c  is  the  sound  speed.  Note 
that  if 


(3.4) 


f  _  J_y  2()^l)qY 

{  M?j  (yfl)u? 


then  (3.1)  has  no  solution.  The  two  values  of  M.2  where  equality  holds  in  (3.4)  correspond  to 
exactly  one  burned  state  each.  These  burned  states  are  the  Chapman  Jouguet  points;  see  Fig.  1. 

For  values  of  M.2  greater  than  the  Chapman-Jouguet  detonation  value,  there  are  two 
possible  burned  states.  The  flow  at  the  unbumed  state  is  supersonic,  as  is  the  case  for  the 
downstream  side  of  a  non-reacting  shock  wave.  At  the  strong  detonation  state  the  flow  is 
subsonic,  which  is  also  normal  for  a  non-reacting  shock  wave.  However,  at  the  weak  detonation 
state  the  flow  is  supersonic,  and  this  violates  the  non-reacting  shock  wave  entropy  condition. 

The  effect  of  the  Mach  number  on  the  flow  pattern  for  (2.7)  at  a  burned  state,  is  as  follows. 
The  linearization  of  (2.7)  at  a  burned  state  has  the  following  eigenvalues: 


(3.5) 


a.  CTi  =  1  - 

b.  02=  1 

C.  O3 


R  RH^  ^  2(u^  -  YP+/P^.) 

m^V^Cp  27U+ 


.4  =  5 1±-\/ '  - 


4Dk0(T*) 


Thus,  the  strong  detonation  state  has  a  two  dimensional  stable  manifold,  and  the  weak  detonation 
state  has  a  one  dimensional  stable  nunfold.  From  this  simple  fact  it  is  already  clear  that  although  a 
weak  detonation  may  be  possible,  .i  strong  detonation  is  much  more  likely. 

The  eigenvectors  correspond.ng  to  O] .  02.  03,  O4  are  as  follows. 
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X,  =  (1,  0.  0,  0) 

X2  =  ((Y-1).  m^v+(l-a,  ).0. 0). 

Y  _  f  n  rr  1 

l^Ym2v+(03  -  <Ji)  L  -  03 


lYin2v,.(04  -  Oi) 


L  -  O3 
L(04  -  1) 
L  -  04 


From  (2.7a)  wc  note  that  Hy  >  0  if  H  -  H4.  +  qe  >  0.  On  the  surface  H  -  H+  +  qe  =  0 
we  have  that 

(3.6)  (H  -  H+  +  qe)y  =  --^r  >  0. 

m-^Cp 

Therefore  the  region  where  Hy  >  0  is  negatively  invariant,  that  is,  the  flow  for  (2.8)  can  only  exit 
this  region;  it  cannot  enter.  On  the  surface  Vy  =  0  (or  Uy  =  0)  we  have 

(3.7) 

Thus,  inside  the  region  Hy  >  0,  the  region  Vy  <  0  is  negatively  invariant.  Also  the  region 
H  ^  H+  is  negatively  invariant  within  Hy  ^  0. 

All  of  the  solutions  that  we  find  will  lie  in  the  region  Hy  ^  0,  H  ^  H.*..  However 
detonation  waves  which  are  close  to  ZND  waves  will  exit  the  region  Vy  ^  0,  attain  a  minimum 
value  of  V  (or  u),  and  maximum  values  of  p  and  p,  and  then  tend  to  the  strong  detonation  state. 
In  order  to  prove  the  existence  of  these  waves,  we  need  another  boundary  of  the  form  V  = 
constant.  A  natural  choice  is 


V  =  Vn  = 


(Y-  1)H. 
(Y+  l)m2v. 


Vq  is  the  value  of  V  at  the  end  state  of  a  non  reacting  shock  wave  with  initial  state  (V. ,  H. );  this 
is  the  minimum  value  of  V  for  a  ZNX)  detonation.  The  part  of  this  surface  lying  within  Hy  >  0, 
e  <  £.,  also  lies  within  Vy  >  0.  Therefore  the  region  V  >  Vg  is  positively  invariant  within  the 
region  Hy  ^  0,  e  <  e. . 

Let  (V,  ,  T,)  be  the  value  of  ( V,  T)  at  the  weak  detonation  state,  and  let  (V*,  T*)  be  the 
value  at  the  strong  detonation  state  Consider  the  region  W  defined  by  Hy  ^  0,  e  ^  e.,  H  ^  H^.. 
V,  ^  V  ^  V.  (See  Fig.  3).  The  flow  for  (2.8)  enters  W  through  a  connected  part  of  the 


& 


.v.v.v, 


boundary,  namely  the  union  of  the  surfaces  V  =  V*,e  =  e, ,  V  =  V.  .  The  flow  leaves  W 
through  two  components,  Hy  =  0  ,  and  H  =  H+  ,  which  are  separated  by  the  positively  invariant 
set  P  defined  by  Hy  =  H  -  H+  =  0,  V*  <  V  ^  Vq.  The  following  theorem  of  Wazewski  [C, 
Wzl,  Wz2],  which  we  quote  from  [C],  implies  that  the  set  of  all  points  in  W  which  eventually 
flow  out  of  W  must  also  have  two  components. 

H 


Fig.  3 


DEFINITION  1.  Let  F  be  a  topological  space  and  let  R  denote  the  real  numbers.  Let  a 
continuous  function  from  F  x  R  F  be  denoted  by  (y,  t)  — >  y*L  This  function  is  called  a  flow 
on  F  if  the  following  conditions  arc  satisfied  for  ail  ye  F  and  s,  t  e  R: 

(a)  y  •  0  =  y 

(b)  y  •  (s  +  t)  =  (y  •  s)  •  t. 

If  F  r)  r  and  R  3  U,  let  F'«  U  be  the  set  of  points  y  •  t  such  that  y  e  T  and  t  e  U. 

DEFINITION  2.  If  F  ^  W,  let  W’  be  the  set  of  points  y  e  W  such  that,  for  some  positive  t, 
y  •  t  e  W.  Let  W~  be  the  set  of  points  y  e  W  such  that  for  any  positive  t,  W  ^  (y  •  (0,  t)).  The 
set  W~  is  contained  in  W  and  is  called  the  exit  set  of  W. 


The  set  W  is  called  a  Wazewski  set  if  the  following  conditions  are  satisfied: 

(a)  Ifye  W  and  cI(W)  d  (y  •  [0,  tj)  then  W  (y*[0,  t]), 

(b)  W"  is  closed  relative  to  W*. 


THEOREM  (Wazewski).  If  W  is  a  Wazewski  set  then  W~  is  a  strong  deformation  retract  of 
W  *  and  W  *  is  open  relative  to  W  *. 

The  set  W  that  we  have  described  above  is  a  Wazewski  set,  as  is  the  union  of  W'  with 
the  set  Q  described  in  the  next  section.  Since  W*  is  homotopic  to  W~,  which  has  two 
components,  W’  must  also  have  two  components. 

What  separates  these  two  components  is  S,  the  stable  manifold  for  the  strong  detonation 
state,  because  any  point  wliich  is  not  in  one  of  these  components  must  stay  in  W  as  y  tends  to 
infinity.  Because  H  is  always  increasing  within  W,  this  solution  must  tend  to  a  rest  state,  namely 
the  strong  detonation  state.  One  boundary  of  this  stable  manifold  is  P;  another  is  F,  the  stable 
manifold  for  the  weak  detonation  state. 

If  there  is  a  connected  set  of  points  in  W  which  tend  to  an  unbumed  state  U,  as  y  tends 
to  -oo ,  and  which  intersects  both  components  of  the  exit  set  of  W,  then  at  least  one  of  these  points 
must  be  in  S.  The  orbit  of  such  a  point  is  a  viscous  structure  profile  for  a  strong  detonation 
wave. 

In  the  next  section  we  give  sufficient  conditions  for  the  existence  of  such  a  set  of  points. 

4.  Some  simple  estimates. 

Orbits  leaving  an  unbumed  state  U  which  lies  in  T  <  Tj  will  stay  in  the  plane  e  =  e. 
until  the  surface  T  =  Tj  is  reached.  Consider  the  curve  G  given  by  Hy  ^  0,  T  =  Tj,  e  =  e.,  and 
Vy  ^  0.  All  of  the  points  of  G  flow  towards  U  as  y  tends  to  -«>.  The  endpoints  of  G  flow 
out  of  Hy  ^  0,  Vy  ^  0,  immediately  .  Our  strategy  is  to  add  a  tunnel  from  W  to  G  so  that  the 
exit  set  of  the  extended  set  is  still  disconnected.  The  existence  of  such  a  tunnel  will  imply  the 
existence  of  a  viscous  structure  profile  from  U  to  the  strong  detonation  state.  The  tunnel  is: 

Q  =  I  (V,  H,  e):  V  ^  V,,  Hy  S  0,  Vy  ^  0,  e.  ^  E  ^  g(T),  and  T  ^  Ti } 

The  function  g  will  be  chosen  so  that  the  flow  enters  Q  through  the  boundary  e  =  g(T),  and  so 
that  g(T)  ^  Ej  for  Tj  ^  T  ^  T*,  where 

There  exist  H  and  V  such  that  1 

Ej  =  minj  £  ^  ^  ^  m2v2/2)/Cp)  =  0  J  ' 


(4.1) 


q(e.-ei)  =  (Cp  -  R/2)(Ti  -  T.) 


See  Fig.  4.  Note  that  e.  -  e,  >  0  whenever  Tj  <  T..  Also  note  that  the  flow  enters  Q  through 
e  =  e.  and  T  =  Tj,  while  it  exits  through  Hy  =  0  and  Vy  =  0.  The  two  components  of  the  exit 
set  are  separated  by  the  point  Hy  =  Vy  =  0,  T  =  T|.  In  order  to  choose  the  function  g,  we  note 
that  within  Q,  we  have  that 

(4  2)  d£  £v  ->-pK(cp(T.  -  T)  -t-  m2(V2 . 

^  '^y  qm2(cp(T  -  T*)  -  m2(V  -  v2)/2  +  qe  -  R(T  -  T*V/V*))  ‘ 

The  right  hand  side  of  (4.2)  is  monotone  in  V  within  Q.  To  see  that  the  denominator  is  monotone 
in  V  for  V*  <  V  ^  V.,  replace  (V*,  H+)  by  (V.,  T.)  and  differentiate  with  respect  to  V,  holding 


Fig.  4.  The  plane  =  H  -  +  qe  =  0. 


T  constant.  This  denominator  is  also  monotone  in  H  as  a  function  of  (H,  V).  In  the  numerator, 
we  decrease  V  to  V,.  In  the  denominator,  we  decrease  V,  holding  T  constant,  until  we  reach  a 
value  Vq  which  satisfies  either  H  =  H+ -  qe  or  Vq  =  V*.  In  the  first  case  we  obtain: 

^  ^  _ XKcp(T>  -  T)4>(T) _ 

dT  -m2qV.(m2(V"5  -VoV*)  -  R(T  -  T,Vo/V*)) 


(4.3) 


^  ^  Xkc^(T,  -  T)4>(T) 
dT  Rin2q2V,(e-f(Vo))’ 


In  the  second  case,  we  decrease  the  denominator  further  by  decreasing  H  while  holding 
V  =  V*,  until  we  reach  H  =  H+  -  qe.  This  again  yields  (4.3),  with  Vq  =  V*  and  a  smaller  value 
of  T  in  the  denominator,  however  this  value  of  T  equals  (H+  -  qe  -  m2v*^/2)/Cp.  We  still  obtain 
(4.4),  with  Vq  =  V*_ 

Suppose  g(T)  ^  Ej  for  Ti<T^T*.  Then  e  within  Q.  Furthermore  f(Vo)^Ei 
within  Q  (see  Fig.  4).  Let 


5^kc2  (T.-T)O(T)  _ 

(4.5)  g’(T)  =  — - >  g(Ti)  =  e. . 

Rm  q  V,  g(T)  -  Ej 

Then  the  flow  enters  Q  through  e  =  g(T)  as  long  as  g(T)  ^  E,.  If  we  solve  (4.4),  we  find  that 
along  E  =  g(T)  we  have: 

T 

/  A  -  \2  _  _  2)^  r  \  l-Z'T'  rr\  X /'T»\  JT’ 


(£  -  Ei)2  =  (E-  -  Ei)2  - 


(y-  I)m2q2v, 


Xk(T,-  T)  (|)(T)  dT  . 


Since  we  need  to  have  e  >  Ej  for  Tj  ^  T  ^  T„  we  require 


(4.7) 


Thus,  (4.7)  is  a  sufficient  condition  for  the  existence  of  a  strong  detonation  structure  profile. 
Since  both  X  and  k  can  depend  on  T  [Wi],  it  is  important  to  keep  them  inside  the  integral. 

For  a  CJ  detonation  wave,  (V,.  T.)  =  (V*,  T*).  Condition  (4.7)  still  implies  the  existence 
of  a  structure  profile. 

As  X  tends  to  zero,  with  the  other  parameters  held  constant,  condition  (4.7)  must  be 
satisfied.  In  particular,  since  e.  >  e,  whenever  T.  <  Tj,  we  see  that  the  strong  and  CJ  detonations 
always  have  viscous  structure  w  hen  X  is  sufficiently  small.  This  is  the  limit  of  the  "small  rate 
parameter"  considered  by  Wood  ( W  o  1 1 

The  region  (Vy  <  0)n(H^  >  V  >  V,)n(£  <  g(T))  is  negatively  invariant,  and  contains 

the  weak  detonation  point,  and  an  interval  of  unbumed  states,  including  one  corresponding  to 
E  =  Ej,  T  =  Tp  in  its  boundary.  Since  the  exit  set  of  this  region  is  homotopic  to  a  circle,  and  the 
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region  itself  is  contractible,  another  application  of  Wazewski's  Theorem  shows  that  the  one 
dimensional  stable  manifold  for  the  weak  detonation  point  must  be  trapped  in  this  region. 
Consequently  one  unbumed  state  with  e  =  e^,  T  =  T^,  V  =  V^,  which  does  not  satisfy  (4.7) 
must  be  connected  to  the  weak  detonation.  Since  this  weak  detonation  is  a  boundary  for  the  stable 
manifold  of  the  strong  detonation  point,  all  unbumed  states  with  £-  >  must  be  connected  to  the 
strong  detonation  point.  As  X.  tends  to  zero,  (e^,  T^,  must  tend  to  (ei,Tj,Vi). 

5.  A  Necessary  Condition 

We  obtain  a  necessary  condition  for  the  existence  of  viscous  structure  using  energy 
estimates  similar  to  those  used  in  work  on  premixed  laminar  flames  [BNS,  Ma,  Wg].  Observe  that 


(5.1) 


H  =  H  -  ^(H^-H)<t)(T). 

’  nrCp 


If  we  multiply  (5.1)  by  1,  H,  and  then  Hy,  and  integrate  each  equation  from  -oo  to  +«»,  we 
obtain,  using  the  fact  that  all  derivatives  tend  to  zero  as  y  tends  to  «»: 


(5.2)  a. 


H+-H.  =  qe.  =  f  i^(H^.H)(t)(T)dy  , 

J 

•  oe 

oo 

oo 

b.  J(Hy)2dy=  f-^(H^-H)H0(T)dy 

-OO  J  rn  Cp 


c. 


OO 

J  (Hy)2  dy  =  f  ^(H^  -  H)<D(T)Hy  dy 


^  Ik 


- (H^  -  H)  <>  ((H  -  m^V  2/2)/Cp)  dH  . 

m-^CpV 


1 


Here  we  have  used  the  fact  that  V.  is  the  maximum  value  of  V  on  any  detonation  structure 
profile.  Ck>mbining  (5.2)(a,  b,  c)  we  find  that; 

(5.3)  (R,  -  H  )2m2cpV.  =  (qY  mrcpV  >  J  2Xk(H^  -  H)0({H  -  m2v2/2)/Cp)  dH  . 
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•  •  • 
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T.  +  qY./Cp 

2XkCp(qY.  -  Cp(T  -  T.))(|>(T)  dT  . 

Wc  may  obtain  a  clearer  interpretation  of  (5.3)  with  a  little  rearrangement  and  change  of  variables: 


(5.3') 


m^V  S 


1 

r  2X)c 

— -<1  -x)<KT.-HqY.T/Cp)  dx 
•0  ^ 


Thus,  (5.3')  constitutes  a  necessary  condition  for  the  existence  of  a  detonation  structure  profile. 
Note  that  if  we  replace  V,  by  V+,  then  we  also  obtain  a  necessary  condition  for  the  existence  of  a 
deflagration  structure  profile.  Also,  (5.3')  is  satisfied  whenever  m^V.  ^  Xk  (J>(T.  +  qYyCp)/Cp. 


6.  Arbitrary  Lewis  Numbers 


When  the  Lewis  number,  L  =  VpDCp  ,  is  not  identically  one,  then  we  must  work  with  all 
four  equations  of  (2.7).  Consequently  the  region  R  must  be  extended  to  a  region  in  R^,  We 
require,  therefore,  additional  boundaries  for  R,  that  is,  upper  and  lower  bounds  for  Y,  It  is 
interesting  that  we  can  obtain  natural  bounds  for  Y  in  terms  of  H,  similar  to  the  bounds  that  have 
been  obtained  for  premixed  laminar  flames. 


Lemma.  Let 


(6.1) 


A*  =  inf(L,  1) 
A*  =  sup(L,  1) 


Then 

(6.2)  A,(H^  -  H)  <  qY  ^  inf(qe,  A*(H+  -  H)) 

defines  a  negatively  invariant  region  for  (2,7). 

Proof.  We  proceed  exactly  as  in  ( Wgl,  using  (2.8): 
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H  -  H+  +  =  (H  -  H+)v  +  -3-Yv 


3^  n-n.^  ~ 


^  (H-H^)y  + 

A* 


=  H-H+  +  qY 

^  H-H++-^Y. 

A, 

This  last  quantity  is  zero  on  the  boundary  of  the  region  defined  by 

H-fT -•--^Y  >  0. 

A* 

Consequendy  this  region  is  negatively  invariant  for  (2.7).  Similarly,  the  region  defined  by 

H-H++-3-Y  S  0 

A* 


is  negatively  invariant,  within  the  region  Y  ^  e.  This  last  condition  is  required  to  ensure  that 
Yy  ^  0.  The  region  Y  ^  e  is  negatively  invariant  because  on  Y  =  e, 

^<e-Y)  =  e,-Yy 

=  Ey  ^  0.// 

The  topological  argument  is  a  little  more  sophisticated  than  the  one  given  in  §3.  The 
region  W  is  defined  as  before,  with  the  additional  inequality  (6.2).  The  exit  set  for  W  is  now 
connected;  using  the  lemma  we  see  that  it  is  a  union  of  four  parts:  Z^,ut  =  Zj  uL2u2^uZ4,  where 

Zi  =  {  (V, H, e,  Y):  H  =  H^.,O^Y<e,0<e<E,Vo<V<V*  } 

L2  =  {  (V,  H,  E,  Y):  Y  =  E,  0  <  E  <  E.,  Vo  <  V  <  V„  0  ^  H+  -  H  S  qE  } 

13  =  {  (V,  H,  E,  Y):  H.,-  H  =  qE,  0  <  E  <  E-,  Vq  <  V  <  V*,  0  ^  Y  ^  E  } 

14  =  {  (V,  H,  E,  Y):  Y  =  0,  0  <  E  <  E.,  Vq  <  V  <  V„  0  ^  H+  -  H  ^  qE  }. 

Note,  however,  that  is  homotopic  to  a  circle.  The  circle  may  be  visualized  as  a  path  from 

Zj  to  Z2  to  Z3  to  Z4  and  back  to  Zj.  This  path  cannot  be  contracted  to  a  point  within  Zq^ 

b«:ause  Z^ut  docs  not  contain  the  positively  invariant  line  segment 


1!^ 


VSf 


P  =  {  (V,H.e.  Y):  Y  =  H-H^  =  e  =  0,  Vo^V^V*}  . 

Again,  Wazewski's  Principle  implies  that  the  set  of  all  points  in  W  which  eventually  flow 
out  of  W  must  also  be  homotopic  to  a  circle.  The  curve  G  of  §4,  given  by  Hy  >  0,  T  =  Tj, 
£  =  £-.  Vy  ^  0,  is  now  extended,  via  (6.2),  to  a  set  G*  which  is  homeomorphic  to  a  disk.  The 
tunnel  Q  is  also  extended  via  (6.2).  The  estimates  of  §4  now  show  that  the  exit  set  of  WuQ  is 
also  homotopic  to  a  circle.  The  boundary  of  G*  is  homeomorphic  to  a  circle,  and  intersects  the 
exit  set  of  WuQ  with  non-zero  degree.  Since  a  disk  is  not  homotopic  to  a  circle,  at  least  one  of 
the  points  of  G*  must  not  exit  W.  This  point  must  in  fact  tend  to  the  strong  detonation  burned 
state.  The  sufficient  condition  (4.7)  is  now 


£.  >  Cj 
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(7-  Um^q^V, 


A.kA*(T,  -  T)<I)(T)  dT 


7.  Behavior 

We  have  shown  that  strong  or  CJ  detonation  structure  profiles  exist  if  (4.7)  is  satisfied. 
From  this  we  can  sec  that  the  existence  of  these  profiles  depends  on  the  values  of  X,  <|),  and  k, 
relative  to  m,  q,  and  7,  between  the  unbumed  state  and  the  weak  detonation  point.  We  will  now 
show  that  the  behavior  of  the  solution  depends  strongly  on  the  values  of  these  parameters  near  the 
strong  detonation  point 

At  the  strong  detonation  point  the  linearization  of  (2.8)  has  two  negative  eigenvalues; 


4Dk({i(T*) 


The  relative  sizes  of  aj  and  03  determine  the  node  structure  of  the  flow  for  (2.8)  at  the  strong 
detonation  point.  If  Dk(})(T*)/u*^  is  very  small,  or  if  L  is  very  small  (which  would  be  physically 
unusual),  so  that  0  >  03  >  Oj  ,  then  the  flow  has  a  node  tangent  to  the  eigenvector 


a. 

Gj  =  -1 

7 

b. 

L 

03  =  j 

q(7  -  1) 
7m^V*(03  -  o,) 


•,  q.  03  -  1 


L(a3  -  1)  ~ 

L  -  03  J 


1 


'  "j* 


Thus,  all  but  one  of  the  structure  profiles  approaches  the  strong  detonation  state  tangent  to  X3. 
The  one  orbit  that  approaches  tangent  to  Xj  is  the  purely  non-reacting  shock  profile  which 
connects  the  weak  detonation  state  to  the  strong  detonation  state.  The  (V,  H,  e)  components  of 
X3  have  the  signs  (+,  +,  -).  Thus,  at  the  end  of  the  wave,  V  and  H  are  increasing  and  e,  which 
must  be  monotone,  is  decreasing.  Since  V*  <  V.,  V  cannot  be  monotone  throughout  the  entire 
solution.  Once  the  solution  leaves  the  region  Vy  <  0,  it  cannot  re-enter,  because  H  is  monotone 
increasing.  Hence  V  attains  its  minimum  value  at  a  single  point  of  the  solution.  This  minimum 
value  must  be  greater  than  Wq,  the  minimum  value  for  a  ZND  detonation. 

The  pressure  must  also  attain  an  extremum  in  this  case; 

dy  V  dy  '  dy 

=  -S-fSf-  m2fn.TiJi.vJl . 


Since  the  flow  is  tangent  to  X3  at  the  strong  detonation  state,  we  have,  as  the  solution  approaches 
that  point, 

=  “Vi .  <r-  1 

Y(O3-0i)  , 

- - V  M-  -2  -  1  -  (r  -  H  M-2, 1 

CpVYVCaj  -  a,)l  H  V  u*2  J  ’) 

CpV*YV*(03- 01)121,  V  u*2  j  j 


-  Y  <0. 


Since  p  is  increasing  when  Vy  <  0,  p  must  attain  a  maximum  in  the  interior  of  the  wave. 
The  temperature  behaves  differently.  Near  the  strong  detonation  state,  we  have 

-  V, 

y  ■  I  J  y 


Thus,  if  M*'^  is  less  than  y.wc  sec  that  T  decreases  near  the  burned  state.  Consequently  T  must 
attain  a  maximwn  value  in  the  interior  of  the  wave.  This  is  consistent  with  the  behavior  of  the 
ZND  wave  corresponding  to  this  case  -  see  [Wi,  p.  194-197]. 

If  M*'^  is  greater  than  y,  then  T  decreases  near  the  burned  state  if  O3-O1  is  sufficiently 
small,  but  positive.  In  the  inviscid  limit  for  this  case,  as  D,  X,  and  p  tend  to  zero,  03  tends  to 
zero  and  T  must  increase  near  the  burned  state.  It  is  a  natural  conjecture  that  T  is  monotone 
throughout  the  wave  in  this  case.  We  say,  therefore,  that  the  flow  is  strongly  subsonic  near  the 
burned  state  if  yM*^  <  1.  If  this  condition  is  satisfied,  then  the  maximum  value  of  the 
temperature  on  the  surface  Vy  =  0  (the  Rayleigh  line)  occurs  at  a  value  of  H  higher  than  for 
on  the  surface  Vy  =  0, 


dT  T,, 

37  =  v(i*y^  )• 


and 


whereas 


dH 

dV 


-^(1  -  m2)  , 
y  -  1 


>0. 


(See  [Wi]).  As  a  consequence,  the  corresponding  ZND  wave,  which  follows  the  intersection  of 
the  surface  Vy  =  0  and  the  plane  Hy  =  0,  must  have  a  monotone  temperature  whenever  the 
burned  state  is  strongly  subsonic. 

One  may  explain  this  phenomenon  by  noting  that  M*  depends  on  the  total  heat  released 
per  unit  mass,  qY.,  and  increases  as  qY.  increases  with  other  parameters,  particularly  m,  held 
constant.  Thus  the  strong  detonation  burned  state  is  strongly  subsonic  if  the  heat  release  is  too 
weak,  relative  to  the  strength  of  the  wave,  to  create  a  temperature  spike.  In  experiments  [see  Wi, 
6.2.1  and  references  cited  therein],  strong  detonation  waves  are  observed  principally  when  a 
piston,  or  other  external  force,  is  used  to  overdrive  the  wave;  in  this  respect  strong  detonation 
waves  resemble  inert  shock  waves.  Detonation  waves  which  are  not  overdriven  will  decay  to  a 
Chapman- Jouguct  detonation.  Chapman  Jouguet  detonations  may  be  thought  of  as  reacting  shock 
waves  that  are  driven  by  the  reaction  with  no  external  force.  Thus  these  waves  are  "pure"  reacting 
shock  waves,  and  qY.  is  a  parameter  between  inert  shock  waves  with  monotone  temperature 
profiles  and  Chapman  Jouguet  waves  with  temperature  spikes. 

As  M*  tends  to  one,  O]  tends  to  zero,  which  leads  us  to  the  next  case.  If  Dk(])(T*)/u*“ 
is  very  large,  or  if  L  is  very  large,  or  if  M*  is  very  close  to  1,  so  that  0  >  Oi  >  03,  then  the  flow 
forms  a  node  tangent  to  the  eigenvector 


X,  =  (1,0,  0,0). 
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In  this  case  the  solution  may  behave  in  an  unusual  manner.  Detonations  with  Y.  close  to  Ej, 
(T.  close  to  T()  will  follow  the  trajectory  of  the  weak  detonation,  and  then  turn  near  the  weak 
detonation  burned  state,  and  approach  the  strong  detonation  point  along  the  trajectory  of  the  inert 
shock  profile.  Thus  these  structure  profiles  look  like  a  weak  detonation  followed  by  an  inert  shock 
wave.  See  Fig.  5,  where  a  heuristic  picture  of  the  flow  in  the  stable  manifold  is  presented^. 
Similar  observations  have  been  made  in  [Wo2,  FD,  LL],  and  most  recently  for  numerical 
calculations  and  for  a  simpler  model  in  [CMR],  where  it  was  noted  that  these  pathological  waves 
are  actually  numerically  stable  as  solutions  of  the  time  dependent  reacting  compressible  Navier 
Stokes  equations.  These  authors  also  observed  (numerically)  an  interesting  phenomenon,  namely 
that  if  the  turning  point  of  these  solutions  is  sufficiently  close  to  the  weak  detonation  state,  that  is, 
if  the  spatial  separation  between  the  weak  detonation  and  the  inert  shock  wave  is  sufficiently  large, 
then  a  bifurcation  occurs  wherein  these  two  waves  decouple  and  the  inert  shock  moves  slower  than 
the  weak  detonation  -  as  predicted  in  [Wo2].  In  this  sense,  a  weak  detonation  can  be  observed. 
Weak  detonations  are  also  observed  experimentally  [FD],  as  a  consequence  of  very  complicated 
chemistry,  change  in  equation  of  state,  and  endothermic  or  reversible  reactions.  A  tixjrc  significant 
point,  made  in  [CMR],  is  that  since  k<{»(T*)/u*^  can  be  significantly  large,  it  is  important,  in 
making  machine  calculations,  to  use  an  approximation  scheme  which  does  not  add  too  much 
artificial  diffusion,  because  this  can  radically  change  the  character  of  the  solutiot . 


(V.H„0) 


Fig.  5.  The  flow  m  the  stable  manifold  when  0  >  Oj  >03. 

^  Note  that  we  have  not  proved  that  the  suble  manifold  has  no  folds,  so  that  V  and  e  are  global  coordinates  for  the 
stable  manifold,  or  that  there  is  a  monotone  rclauonship  between  Y.  and  the  solution  curves,  such  as  depicted  in 
Fig.  5.  For  a  proof  of  such  monotoniciiy  for  premixed  laminar  flames  with  L  >  1  see  [Ma). 
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For  larger  values  of  Y. ,  there  will  be  one  wave  which  approaches  the  strong  detonation 
state  tangent  to  X3  ,  and  others  which  approach  tangent  to  X]  but  with  V  increasing.  For  the 
singular  wave  tangent  to  X3  all  quantities  will  be  monotone,  but  for  the  others  the  pressure  and 
density  will  attain  maximum  values.  Since  H  =  CpT  +  (m/p)^/2  is  constant  along  Xj  ,  the 
existence  of  a  density  peak  in  these  waves  implies  the  existence  of  a  temperature  peak. 

For  the  O  detonation,  the  burned  state  Mach  number  is  1,  so  that  the  above  remarks  apply 
to  the  behavior  of  this  wave.  All  but  one  of  the  Cl  detonation  structure  profiles  must  approach 
tangent  to  the  eigenvector  ±(1,0,  0,  0).  However  solutions  cannot  approach  along  +{1,  0,  0,  0) 
(the  right  side)  because  Vy  >  0  there.  Consequently  there  is  one  monotone  CJ  structure  profile 
(presumably  the  one  with  minimum  heat  release)  and  the  rest  have  peaks  in  pressure,  density,  and 
temperature. 

8.  The  ZND  Limit 

We  have  noted  that  as  X  tends  to  zero,  (or  X  and  D,  with  L  constant),  then  the 
sufficient  condition  (4.7),  or  (6.3),  must  be  satisfied.  The  structure  profiles  satisfy  natural  a  priori 
bounds,  namely 

Vo^V^V., 

H.  ^H^H+, 

O^Y^e^Y.. 

Thus,  for  any  sequence  X.n  0,  there  is  a  corresponding  sequence  (V^,  H^,  Gn*  ^n) 
structure  profiles  with  fixed  end  states,  and  which  are  uniformly  bounded.  Since  H,  e,  and  Y  arc 
monotone,  and  V  has  at  most  one  minimum,  these  profiles  are  also  uniformly  bounded  in  total 
variation.  By  Helly's  theorem,  some  subsequence  of  this  sequence  must  converge  to  a  limit  in 
l/oc  •  Taking  a  further  subsequence  we  obtain  convergence  a.e.  The  limit  function  is  therefore  a 
weak  solution  to  (1.1  a,b,c,d').  If  T.  <  T,  then  the  limit  is  a  ZND  strong  or  CJ  detonation.  If 
T-  =  Ti  then  the  limit  is  a  continuous  weak  detonation  which  follows  the  curve 
(Hy  =  0)n(Vy  =  0)n(e  =  Y)  from  (Vj,  Hj,  Ej,  Yj)  to  (V*,  H„  0,  0). 

9.  The  Second  Law  of  Thermodynamics 

In  several  shock  structure  problems,  and  particularly  in  MHD[CS1,  CS2,  Ge],  the  entropy 
flux  has  played  an  important  role.  In  MHD  the  structure  equations  take  the  form 


.  -  «  -  ■  •  ^  ' 


where  B  is  a  positive  diagonal  matrix  and  P  is  the  entropy  flux.  The  gradient  like  structure  that 
this  gives  to  the  problem  is  essential  to  our  understanding  of  the  solution  to  this  very  complicated 
system.  In  other  areas  it  is  useful  to  postulate  that  the  entropy  flux  must  be  monotone 
[HaSel,  HaSe2];  the  inequality  expressing  this  monotonicity  is  called  the  Clausius -Duhem 
inequality  .  We  show  here,  assuming  only  that  the  reaction  is  exothermic  and  irreversible 
(r(p,  Y,  T)  2:  0),  that  the  entropy  flux  is  monotone  along  a  viscous  structure  profile  for  a  plane 
steady  detonation  or  deflagration  wave.  We  have  not  used  this  fact  in  the  above  discussion, 
although  it  would  have  proved  useful  if  we  did  not  already  know  that  H  must  be  monotone. 

The  entropy  flux  is 

P  =  mS  -  XT^/T. 

Using  Gibb's  law;  TdS  =  de  +  pdV  -  qdY  [Wi],  we  find  that 


dP 

ar 


-  Ef*  dJV  _  d  Y'l  ^ 

■  TVdx  P  dx  ‘laxj^  t2  T 


Using  (1.2),  we  have 
dP  XtI  1 


ar 


=  -t([[p(u^/2  +  e(p.  T,  Y))  +  p(p,  T)]u]^  -  (puuj,  -  (qpDY;^);^) 


^mfde^  dV  dY 


XTi  pul 


dP 


Note  that  only  vanishes  at  rest  points  of  (1.2).  However,  (1.2)  is  not  gradient  like  with 
respect  to  P,  because 

dP  =-^dT  +  2^  -^du; 


since  q/T  does  not  vanish,  neither  does  the  gradient  of  P.  This  noay  be  an  anomaly  which  is  due  to 
the  simple  representations  of  the  chemistry  and  the  reaction  rate  which  are  used  here. 
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